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Abstract 

In the framework of quantum field theory we discuss the emergence of a phase 
locking among the electromagnetic modes and the matter components on an 
extended space-time region. We discuss the formation of extended domains 
exhibiting in their fundamental states non-vanishing order parameters, whose 
existence is not included in the Lagrangian. Our discussion is motivated by the 
interest in the study of the general problem of the stability of mesoscopic and 
macroscopic complex systems arising from fluctuating quantum components 
in connection with the problem of defect formation during the process of non- 
equilibrium symmetry breaking phase transitions characterized by an order 
parameter. 

P.A.C.S.: ll.lO.-z, 64.90.+b, ll.30.Qc 



I. INTRODUCTION 

Complex systems made up of quantum components are usually remarkably stable at 
mesoscopic and macroscopic space-time scales. On the other hand, quantum fluctuations are 
the dominant feature at the microscopic scale of the quantum components. The necessity of 
taking into account such a double feature is reflected in the usual quantum field theory (QFT) 
prescription that the Lagrangian of the complex system built upon the quantum fields should 
be invariant under the local phase transformation of the quantum component field ^/'(x, t) — >■ 
ip'l'x,!) = exp{ig6{:s.,t))ilj{'x,t). Local phase invariance is the QFT solution to the problem 
of building a stable system out of fluctuating components. The requirement of local phase 
invariance demands the introduction of gauge fields, e.g. the electromagnetic (e.m.) field 
A^(x, t), such that the Lagrangian be also invariant under the local gauge transformation 
yl^(x, t) A^(x, t) — 9^6'(x, t). Such a transformation is devised to compensate terms 
proportional to dfj_6{'x,t) arising in the Lagrangian from the kinetic term for the matter 



field ilj{:K,t). This is a well known story. In the present paper, given the above connection 
between the matter field and the e.m. field, we wish to discuss, in the frame of QFT, 
the role played by the e.m. field in the locking of the phases of the e.m. modes and of 
the matter components on an extended space-time region. Furthermore, we will discuss 
extended domains exhibiting in their fundamental states non-vanishing order parameters, 
whose existence is not included in the Lagrangian. 

The interest in the general problem of the stability of mesoscopic and macroscopic com- 
plex systems arising from fiuctuating quantum components also finds one strong motivation 
in the study of the physically relevant problem of defect formation during the process of non- 
equilibrium symmetry breaking phase transitions characterized by an order parameter^. A 
topological defect may indeed appear in such a process whenever a region, surrounded by 
ordered domains, remains trapped in the "normal" or symmetric state. Examples of topo- 
logical defects are vortices in superconductors and superfluids, magnetic domain walls in 
ferromagnets, and many other extended objects in condensed matter physics. On the other 
hand, topological defects, such as cosmic strings in cosmology, may have been also play- 
ing a role in the phase transition processes in the early Universe^. The phenomenological 
understanding of the defect formation in phase transitions is provided by the Kibble-Zurek 
scenario^'^. By considering the surprising analogy between defect formation in solid state 
physics and in high energy physics and cosmology^, it has been also stressed that the anal- 
ysis of the formation of defects in phase transitions becomes a "diagnostic tool"^ in the 
study of non-equilibrium symmetry breaking processes in a wide range of energy scales. 
Questions such as why extended objects with topological singularities are observed only in 
systems showing some sort of ordered patterns, why defect formation is observed during 
the processes of phase transitions, why the features of the defect formation are shared by 
quite different systems, from condensed matter to cosmology, etc., have been specifically 
addressed in refs.^'^ and the dynamics of defect formation has been extensively studied in a 
large body of literature in QFT; see, as general refs,^ and^. 

In these studies, in dealing with the presence of a gauge field in the process of spontaneous 
symmetry breakdown a crucial role is played by the well known Anderson-Higgs-Kibble 
(AHK) mechanism^*^'^^ , where the gauge field is expelled out of the ordered domains and 
confined, through self- focusing propagation, into "normal" regions, such as the vortex core, 
having a vanishing order parameter, i.e. where the long range correlation modes (the Nambu- 
Goldstone modes) responsible for the ordering are damped away. In the present paper, going 
beyond the well established AHK mechanism, our attention is focused on the dynamics 
governing the radiative gauge field and, as said above, its role in the onset of phase locking 
among the e.m. modes and the matter components. In the AHK mechanism the gauge 
field removes the order in the regions where it penetrates, thus describing the self-focusing 
gauge field propagation in ordered condensed matter as well as in asymmetric vacuum in 
elementary particle physics. Here we study the role of radiative gauge field in sustaining the 
phase locking in the coherent regime. 

We choose as our model system an ensemble of a given number N of two-level atoms, 
which may represent rigid rotators endowed with an electric dipole. The interaction of 
these atoms with the e.m. quantum radiative modes will be considered. Moreover, we will 
examine the effects on the system of an electric field generated by an external source or else 
by an impurity introduced in the system, thus making contact with the family of the so 
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called Jaynes-Cummings models^ , extensively studied in the literature in connection with 
quantum optics problems (see e.g. ref.^'^ for detailed analysis). Our discussion in the present 
paper may be indeed of some relevance to quantum optics, as well. As a matter of fact, a 
system of N two-level atoms interacting with quantized e.m. modes has been known since 
long time to provide from a formal point of view, under some convenient approximations 
and restrictions, a strong analogy between the laser phase transition and the conventional 
phase transitions in spontaneously broken symmetry theories^^'^^, although the meaning of 
the constants in the potential function for the order parameter is different. The key point in 
such a connection between coherent (laser) light, the N atom system and phase transition 
is in the observation that, under convenient conditions, the behavior of the e.m. mode is 
described by the potential 

V{u, u*) = -a\u\' + p\u\' + ^ ^ p{\u\' - , (1.1) 

where, in the Haken notation (see eq. (2.19) of^"^ or eq. (VI. 4. 24) of^^) a and P, with P > 0, 
are convenient coefficients and u denotes the classical e.m. amplitude corresponding to the 
quantum e.m. field amplitude (in the interaction representation the electric field strength is 
decomposed as E — uex'p(—iu}ot)+u* exp{iuJot) where u is the slowly varying amplitude and 
u!o the atomic resonant frequency). The essential point is that the (mean value of the) "order 
parameter" u minimizing the potential V{u,u*) is zero (disordered or symmetric state) for 
a < and non-zero for o; > 0, with \u\^ = ^ (ordered or asymmetric state). In this 
latter case the system is said to be above threshold (the threshold is set at ck = 0), i.e. it is 
lasering. Of course, in the Lagrangian formalism the coefficient (—a) denotes the "squared 
mass" of the field, whose sign, as well known, controls the occurrence or not of spontaneous 
symmetry breakdown. In the Haken interpretation a is the pump parameter whose tuning 
may carry the system far from the equilibrium, i.e. in the lasering region. Thus in the phase 
transition between the disordered and the ordered state the order parameter u changes in 
time from zero to a value proportional to y/a. In his analysis, Haken also considers the 
Hamiltonian in the interaction representation 

H = n-f{b^S- + bS+) , (1.2) 

which is a Jaynes-Cummings-like Hamiltonian, indeed. In eq. (1.2) 7 is a coupling constant 
which is proportional to the atomic dipole moment matrix element and to the inverse of 
the volume square root V~^^'^, b is the e.m. quantum field operator (associated to the c- 
number amplitude u), are the atomic polarization operators. In the Haken discussion 
the atomic variables are integrated out at some point of the computation since his interest is 
mostly focused on the e.m. lasering effect. In our following analysis, instead, we keep them 
and show that the phase locking between them and the e.m. mode can be reached under 
convenient boundary conditions. 

II. THE MODEL 

Let us start by assuming that transitions between the atomic levels are radiative dipole 
transitions. We thus disregard the static dipole-dipole interaction. Moreover, the system 
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is assumed to be in a thermal bath kept at a non-vanishing temperature T. Under such 
conditions the system is invariant under dipole rotations. We use natural units h = 1 = c. 
We assume the system be spatially homogeneous and denote by N the number of atoms per 
unit volume. The N atom system may be collectively described by the complex dipole wave 
field 0(x, t). In Section IV we will also use the known formal equivalence (see e.g. Section 
III. 6 of ref.^'') of the system of two-level atoms with a system of | spins. The dipole wave 
field 0(x, t) integrated over the sphere of unit radius r gives: 

I dn\<j>i^,t)\' = N , (2.1) 

where dfl — smOdOdcf) is the element of solid angle and (r, 9, (f) are the polar coordinates of 
r. By introducing the rescaled field x(x, = -^(f)(yi,t) eq. (2.1) becomes 

Jdn\x{^,t)f = l. (2.2) 

Since the atom density is assumed to be spatially uniform, the only relevant variables 
are the angular ones. Thus, in full generality, we may expand the field x(x, t) in the unit 
sphere in terms of spherical harmonics: 

xM^J2^Ut)Yrie,<i>) , (2.3) 

l,m 

which, by setting a^rn{t) = for / 7^ 0, 1, reduces to the expansion in the four levels (/, m) = 
(0, 0) and (1, m),m = 0, ±1. The populations of these levels are given by N\ai^rn{t)\'^ a-nd at 
thermal equilibrium, in the absence of interaction, they follow the Boltzmann distribution. 
Moreover, the dipole rotational invariance implies that there is no preferred direction in the 
dipole orientation, which means that the amplitude of «i,m(i) does not depend on m, and 
that no permanent polarization may develop for such a system in such conditions, i.e. the 
time average of the polarization along any direction n must vanish. We thus write 

ao,o(t) = ao{t) = Ao{t) (^'^^ , 

a^,^{t) = A^{t) e^^i-W e— °* = a^,m{t) e"^""* , (2.4) 

where ai^m{t) = ^i(0 e^s-i,rn{t) _ AQ{t), A-i{t), 5o(t) and 5-i^jn{t) are real quantities. In eqs. 
(2.4) we have also used ujq = where / denotes the moment of inertia of the atom, which 
gives a relevant scale for the system: ujq = k — (note that the eigenvalue of |y on the 
state (l,m), being the squared angular momentum operator, is ^^^^P ~ 1 ~ ^o)- By 
setting the z axis parallel to n and using the explicit expressions for the spherical harmonics 




sin e e''!' ^ - [Ff , (2.5) 
we find 
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= -^^oW^iW cos(c^ - uJo)t , (2.6) 

where cut = 5i,o(^) — ^oit) and whose time average is zero, as it should be. This confirms 
that the three levels (l,m), m = 0, ±1 are in the average equally populated under normal 
conditions and that, as said above, we can safely write J2m |Q^i,m(^)P — 3 |ai(i)p. On the 
other hand, the normalization condition (2.2) gives, at any time t, 

l«o,o(t)|' + E l«i,m(t)r = |ao(t)r + 3 |ai(t)|2 = 1 , (2.7) 

m 

By defining Q as 

Q = \ao{t)\^ + 3 \ai{t)\\ (2.8) 



we thus see from eq. (2.7) that 



|q = . (2.9) 



I.e 



|KWI^ = -i|KWp. (2.10) 

Due to eq. (2.1) (and the rescaling adopted for x(x, t)), eq. (2.9) expresses nothing but the 
conservation of the total number N of atoms; it also means that, as shown in eq. (2.10), due 
to the rotational invariance, the rate of change of the population in each of the levels (1, m), 
m — 0, ±1, equally contributes, in the average, to the rate of change in the population of 
the level (0, 0), at each time t. Consistently with eq. (2.7), in full generality we can set the 
initial conditions at i = as 

|ao(0)|2 = cos2^o , |ai(0)|2 = ^sin^^o , < < | ■ (2.11) 

We exclude the values zero and | since they correspond to the physically unrealistic con- 
ditions for the state (0,0) of being completely filled or completely empty, respectively. By 
properly tuning the parameter 6q in its range of definition one can adequately describe the 
physical initial conditions. For example, = f describes the cquipartition of the field modes 
of energy E{k) among the four levels (0,0) and (l,m), |ao(0)p ~ |ai_^(0)p, m — 0, ±1, as 
typically given by the Boltzmann distribution when the temperature T is high enough, 
ksT 3> E{k). As we will see, however, the lower bound for the parameter Oq is imposed by 
the dynamics in a self-consistent way. 

The possibility of obtaining a non-zero permanent polarization, and thus the dipole 
ordering in the system, is crucially conditioned by the ratio between the populations in the 
atomic levels. Namely, suppose that the atom system is under the influence of an electric 
field E due, e.g., to an impurity, or to any other external agent. Assume E to be parallel to 
the z axis. Then the term 
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(2.12) 



where d is the electric dipole moment of the atom, will be added to the system energy and 
will break the dipole rotational symmetry. It will produce the mixing between the states Yq 
and Y^: Y^ Y^ cos r + Y^ sin r and Y^ -Y^ sin r + Y^ cos r, with 



~m ■ ^ ^ ^ 



Due to the mixing thus induced, the polarization is now given by 



+ ^AQ{t)Ai{t) cos2t cos([ uj - y/uf+m^ ]t) , (2.14) 
V 3 

to be compared with eq. (2.6) and whose time average is non-zero: Pn — -^{Al — Af) sin 2t. 
The non-zero difference in the level populations {A^ — Af) is therefore crucial in obtaining 
the non-zero polarization. We will study under which conditions such an occurrence can be 
realized. 

In the following, we restrict ourselves to the resonant radiative e.m. modes, i.e. those 
for which k — ^ — cuq, and we use the dipole approximation, i.e. we put exp(ik • x) « 1 
in our formulas, since we are interested in the macroscopic behavior of the system. This 
means that the wavelengths of the e.m. modes we consider, of the order of are larger (or 
comparable) than the system linear size. Let Cr(k, t) denote the radiative e.m. field operator 
with polarization r and Ur(k,t) = -^Cr{k,t) the rescaled one. The field equations for our 
system are 



13,16. 



l- 



dt 21 ^) ~ ^ S ^ ^ V 2 



durCk, t) 



K(k,i)e-'=*-i.t(k,i)e^i x(x,t), 

^ dVP \/f e''* / c^O(6.-x)|x(x,t)|2 , (2.15) 



where d is the magnitude of the electric dipole moment, p = y and is the polarization 
vector of the e.m. mode (for which the transversality condition k • = is assumed to 
hold) . Wc remark that the enhancement by the factor ^/N appearing in the coupling d 
in eqs. (2.15) is due to the the rescaling of the fields. We will comment more on this point 
later on. 



III. THE FIELD EQUATIONS AND THE PHASE LOCKING 

By resorting to the discussion of the eqs. (2.15) presented in ref.^^, our task is now to 
analyze the implications of these equations with reference to the role played by the e.m. 
modes in the onset of the phase locking between them and the dipole field. 
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First, let us observe that use of eqs. (2.3) (and (2.4)) in (2.15) gives the set of coupled 
equations 

aoit)^QYl <W (3.1) 

m 

ai,m(t) = -Q Um(t) ao(t) (3.2) 
ii^it) = 2 n a*(t) a,,m{t) , (3.3) 

where ai^rnit) = Oii^mit) e*"^"* (cf. eq. (2.4)), the dot over the symbols denotes the time 
derivative, Um is the amplitude of the e.m. mode coupled to the transition (l,m) ^ (0,0) 

Eqs. (3.1)-(3.3), as well as the eqs. (2.15), from which they have been derived, appear to 
be not invariant under time-dependent phase transformations of the field amplitudes. We 
want to investigate how gauge invariance can be recovered. 

Use of the conjugate of eq. (3.2) in (3.1) gives, consistently with eq. (2.9), the conserva- 
tion law Q — 0. Moreover, use of the conjugate of eq. (3.2) in (3.3) leads to 

^K(i)r = -2||ai,^(i)r (3.4) 

Since the amplitude |Q;i,m(^)| = |cii,m(t)| does not depend on m, eq. (3.4) shows that also the 
amplitude 11*^(^)1 does not depend on m. Eq. (3.4) shows the existence of another constant 
of motion; namely, putting \u{t)\ = \ujn{t)\ and using |ai(t)| = |ai,m(^)|, we can write 

|M(i)|2 + 2 |ai(i)|2 = ^sin^^o , (3.5) 

for any t, where we have also used the initial condition (2.11) and set 

\u{0)f^O. (3.6) 

We observe that since |M(t)P > eq. (3.5) imposes |ai(t)P < | sin^ 6'o and therefore 
|ao(t)P > cos^^o due to (2.7). Note that (3.5) gives (cf. (2.11)) 

\u{t)\' = 2{ |ai(0)r-|ai(t)n , (3.7) 

for any t, which, by exploiting (2.11), reads 

\u{t)\'^^{ \ao{t)\' - cosHo ) . (3.8) 

It is also useful to write 

Um(t) = C/(t)e'<^-W , (3.9) 

with U{t) and (pm{t) real quantities. 

By combining eqs. (2.4) and (3.9) with (3.1)-(3.3) and equating real and imaginary 
parts, we get 
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where we have put 



Ao{t) = nU{t)Ai{t) cosam{t) , (3.10) 

Ai{t) = -nU(t)Ao(t) cos am(t) , (3.11) 

U{t) = 2nAo{t)Ai{t) cosam{t) , (3.12) 

^■Ut) = 2Q^^^^^sina„(t) , (3.13) 



= (^l.mlt) - ^"o(^) - ^mit) ■ (3-14) 



Similarly, we can derive equations for Si^m and So- 

From eqs. (3.10)-(3.12) we see that since their left hand sides are independent of m, 
so the right hand sides have to be, i.e. either cos am{t) = for any m at any t, or is 
independent of m at any t. In both cases, eq. (3.13) shows that (f^ is then independent 
of m, which in turn implies, together with eq. (3.14), that Si^mit) is independent of m. 
Phases thus turn out to be independent of m. We will therefore put (p = (pm, Si{t) = Si^mit) 
and a = We observe that in general the phases can be always changed by arbitrary 
constants. The independence of m of the phases is dictated by the requirement to not violate 
the gauge invariance. Should exist a difference between the phases having different m, this 
difference could be changed by a rotation of the axes and would spoil the gauge invariance. 
In the present case, the independence of m of the phases is of dynamical origin and we will 
find that the phase locking among (5o(t), Si{t) and ip{t), has indeed the meaning of recovering 
the gauge invariance. We will discuss this point in the subsection II. B. 

Summarizing, we can now write u{t) = Ujn{t) and ai{t) = ai^^(t) and from the equations 
(3.1)-(3.3) we get the known-'^''' set of equations: 

do{t) = 3n u*{t) ai{t) (3.15) 
di(t) = -Q u{t) ao{t) (3.16) 
u{t) = 2 Q al{t) ai{t) . (3.17) 

Eqs. (3.15)-(3.17) are fully consistent with the original normalization condition (2.2) (or 
(2.7)), with the conservation (2.9) and with the dipole rotational invariance expressed by 
the zero average polarization (cf. eq. (2.6)). In eq.(3.15) the rate of change of the amplitude 
of the level (0, 0) is shown to depend on the couphng between the levels (1, m), m = 0, ±1 and 
the radiative e.m. mode of corresponding polarization. Each of these couplings contribute in 
equal measure, due to rotational invariance, to the transitions to (0,0). Similarly, in eq.(3.16) 
the rate of change of the amplitude of each level (1, m) is shown to depend on the coupling 
between the the level (0, 0) and the corresponding radiative e.m. mode. Finally, in eq.(3.17) 
the transitions (0,0) <-> (l,m),m = 0,±l control the rate of change of the amplitude of the 
radiative e.m. mode of corresponding polarization. These equations thus reflect the correct 
selection rules in radiative and absorption processes^^^^". Eq. (3.15) describes the fact that 
each of the levels (1, m) may find in the e.m. field the proper mode to couple with, in full 
respect of the selection rules. In this sense, the field concept, as a full collection of e.m. 
modes with all possible polarizations, is crucial here. As already said, eqs. (3.15)-(3.17) are 
fully consistent with the physical boundary conditions and the motion equation (2.15) from 
which they are derived. 
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A. The ground state 



We want to study now the vacuum or ground state of the system for each of the modes 
ao(t), ai{t) and u{t) described by (3.15)-(3.17) (i.e. by (2.15)). 

It is convenient to differentiate once more with respect to time both sides of (3.15). By 
using (3.16) and (3.17) and the constants of motion (2.7) and (3.5) we ehminate the variables 
ai{t) and u{t). We thus find: 

ao{t) = 4 Q^^^Maoit) - 4 Q>o(*)|'ao(i) , (3.18) 
where 70(6*0) = |(1 + cos^ 6*0). Eq. (3.18) can be written in the form 

r 

ao{t)^-YTVMt):al{t)] , (3.19) 

where the potential V[ao{t), aQ{t)] is 

V[ao{t),al{t)] = 2n\\ao{t)\'-j'Mr . (3.20) 

In a standard fashion (see e.g.^^) we may adopt the semiclassical ('mean field') approxima- 
tion in order to study the ground state of the theory. We thus search for the minima of the 
potential V. Let ao,_R(t) and ao,/(t) denote the real and the imaginary component, respec- 
tively, of the ao{t) field: |ao(t)P = ^o(^) — ^o,ii(^) + '^oji^)- The potential has a relative 
maximum at oq = and a (continuum) set of minima given by 

|ao(t)|' = ^(l + cos2^o)=7o(^o) ■ (3.21) 

These minima correspond to the points on the circle of squared radius Jq{0o) in the 
(Qo,i?(^)) Cio,i{t)) plane. Wc thus recognize that we are in the familiar case of a theory where 
the cylindrical 5*0(2) symmetry (the phase symmetry) around an axis orthogonal to the 
plane {ao^ii{t), aoj{t)) is spontaneously broken. The order parameter is given by 70(^0). 
Note that eq. (3.21) does not fix the (ground state expectation) value of the phase field 
6o{t). The points on the circle represent (infinitely many) possible vacua for the system and 
they transform into each other under shifts of the field 60: 6q ^ 6q + a (SO (2) rotations in 
the {ao,R{t),aoj{t)) plane). The phase symmetry is broken when one specific ground state 
is singled out by fixing the value of the Sq field. 

By proceeding as usual in these circumstances^^, we transform to new fields: Ao(t) 
A'(^{t) = Ao{t) — 7o(^^o) and 6^(1) — So{t), so that A'^lt) = in the ground state for which 
Ao{t) = 7o(^^o). Use of these new variables in V leads to recognize that the amplitude Ao(t) 
describes a quasi-periodic mode with pulsation mo = 2^2^' (1 + cos^ ^^q) (a 'massive' mode 

with real mass (1 -|- cos^ ^0)) and that the field 5o(t) corresponds to a zero- frequency 
mode (a massless mode) playing the role of the so-called Nambu-Goldstone (NG) field or 
collective mode implied by the spontaneous breakdown of symmetry. 

We note that when eq. (3.21) holds, use of eqs. (2.7) and (3.5) gives A^ — |sin^^o and 

u'^ = ^sin^^O) moving away from the initial condition values (2.11) and (3.6), respectively. 
In this respect, it is remarkable that the value Qq — 0, which we have excluded in our initial 
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conditions, cf. eq. (2.11), on the basis of physical considerations, consistently appears to 
be the relative maximum for the potential, and therefore an instability point out of which 
the system (spontaneously) runs away. Moreover, as already observed, use of the constant 
of motion laws (2.7) and (3.5) shows that |aoP = implies C/^ = — |cos^^o which is not 
possible since U is real. Finally, we remark that the bound |ao(t)P > cos^ 9 discussed above 
(see the comment after eq. (3.6)) is consistently satisfied by |ao(^)P = 7o(^o) (see eq. (3.21)). 

We consider now the time derivative of both sides of (3.16) and by simple manipulations 
we find the following equation for the amplitude ai{t): 

d^{t) = -a^ai{t) + 12 Q>i(t)|V(^) , (3.22) 

where = 2 f2^(l + sin^^o)- The potential from which the r.h.s. of eq. (3.22) is derivable 
is 

V[a,{t),al{t)]^a'\a,{t)\'-6n\\a^{t)\Y . (3.23) 

In this case there is a relative minimum at ai = and a (continuum) set of relative maxima 
on the circle of squared radius 

\a,{t)\' ^ ^{1 + smHo) = jlieo) . (3.24) 

Note that, for |ai(t)p = 7i(6'o), U"^ = — | cos^ 9q < 0, which is not acceptable since U is real. 
The values on the circle of radius 7i(6'o) are thus forbidden for the amphtude Ai. This is 
consistent with the intrinsic instability of the excited levels (l,m). We have also seen that 
the conservation law (3.5) and the reality condition for U require that |ai(t)P < Isin^^o 
which lies indeed below ^f{9Q), and we note that the value | sin^ taken by when 
|'^o(^)P = 7o(^o) also lies below the bound. In conclusion, the potential V[ai{t),al{t)] 
involved in the dynamics must be lower than | sin^ 9. 

This is enough about the consistency between eqs. (3.18) and (3.22). As men- 
tioned above, we exclude that the amplitude Ai be constantly zero (at the minimum of 
V[ai{t),al{t)]), since this would correspond to the physically unrealistic situation of the 
(0,0) level completely filled. In conclusion, within these dynamical bounds, the field ai{t) 
described by Eq. (3.22) is a massive field with (real) mass (pulsation) a"^ — 2 Q^(l-|-sin^ ^o)- 

Finally, we focus on the e.m. mode u{t) and consider eq. (3.17). By proceeding as above 
by differentiating once more with respect to time we find 

u(t) = -iJ,\{t) - 6 n'^\u{t)\\{t) , (3.25) 

where ji'^ — 2 Vt^ cos29q. The r.h.s. of eq. (3.25) is derivable from the potential 

V[u{t),u*{t)] = i/\u{t)\^ + 2>n^\u{t)\^ + \n'^cos^29Q 

= 3Q'(|^i(i)|' + Jcos2eo)' , (3.26) 

and we note that V [u{t) ^ u* {t)] is nothing but the potential for the e.m. mode given in eq. 
(1.1) for —a = and /3 = 3 Vl?. We are in the case of a theory where the symmetry can 
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be spontaneously broken or not, according to the negative or positive value of the squared 
mass /i^ of the field (the pump in the Haken interpretation), respectively. 

Again, in the semiclassical approximation we search for the minima of the potential 
V[u{t),u*{t)] and see that //^ > for 9o < ^ and the only minimum is at uq — 0. Eq. (3.25) 
then describes quasi-periodic modes with pulsation /i — Q-\/2~cos2^, typically expected 
for a paraboloid potential V[u(t),u*{t)] with cylindrical S0{2) symmetry about an axis 
orthogonal to the plane {u[i{t),uj{t)) and minimum at uq = 0. Here Ufi{t) and Ui{t) denote 
the real and the imaginary component, respectively, of the u{t) field. In such a case we 
have the symmetric solution with zero order parameter uq — 0. This solution describes the 
system when the initial condition, eq. (3.6), holds at any time. This occurrence is, however, 
not consistent with the dynamical evolution of the system moving away from the initial 
conditions exhibited by eq. (3.18), as mentioned above. Luckily, consistency is dynamically 
recovered provided 6'o > f . In such a case, indeed, ii^ = 2 ft^ cos 29o < and the potential 
has a relative maximum at = and a (continuum) set of minima given by 

\u{t)\' = cos 2^0 = = ^'(^o) , Oo>^. (3.27) 

The fact that in the present case = is a maximum for the potential means that the 
system dynamics evolves away from it, consistently with the similar situation noticed above 
for the ao mode where the system spontaneously evolves away from the initial conditions. 
The symmetric solution of the minimum ai uq — is thus excluded for internal consistency 
and the lower bound | for is thus dynamically imposed in a self-consistent way. 

In eq. (3.27) the minima are the points of the circle of squared radius v^(^o) in the 
{uR{t),ui{t)) plane. As in the case of the amplitude analyzed above, the points on the 
circle represent (infinitely many) possible vacua for the system and they transform into each 
other under shifts of the field y?: ^ + a. For > f the phase symmetry is broken, the 
order parameter is given by v{9o) ^ and one specific ground state is singled out by fixing 
the value of the (/? field. 

As usual"*^^, we transform to new fields: U{t) U'{t) = U{t) —v{9q) and v^'(t) ip{t) so 
that in the ground state U'{t) = 0. Use of these new variables in V [u{t) , u* (t)] shows that 

the amplitude U'{t) describes a 'massive' mode with real mass ^J2\n'^\ = 2Q^J\ cos26'o| (a 
quasi-periodic mode) and that the field (p'{t) is a zero- frequency mode (a massless mode). 
This field, also called the "phason" field^^, plays the role of the Nambu-Goldstone (NG) 
collective mode in the theories where symmetry is spontaneously broken. When eq. (3.27) 
holds, it is Af = ^ which lies below the upper bound | sin^ 9 provided 9 > ^. Similarly, eq. 
(3.27) implies = | which satisfies the constraint of being greater than cos^ 6* for 6* > |. 

In conclusion, the e.m. field, as an effect of the spontaneous breakdown of the phase 
symmetry (eq. (3.27)) gets a massive component (the amplitude field), as indeed expected in 
the Anderson-Higgs-Kibble mechanism, and there is also a (surviving) massless component 
(the phase field) playing the role of the NG mode. In the following we show that such a 
massless component is crucially involved in the phase locking of the e.m. and matter fields. 

The emerging picture is then the following. The system may be prepared with initial 
conditions dictated by the conservation of the particle number and given by eqs. (2.11) and 
(3.6), where the value of the parameter 9q is in principle arbitrary within reasonable physical 
conditions. According to the field equations (2.15), the system then evolves towards the 
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mininiTim energy state where |ao(t)P 7^ as in eq. (3.21) and the amplitTide departs 
from its initial zero value. This implies a succession of (quantum) phase transitions^^ from 
the initial uq = symmetric vacuum to the asymmetric vacuum |M(t)P 7^ 0, which means 
that in eq. (3.26) 6*0 has to be greater than |. In this way the lower bound for ^0 is 
dynamically fixed and the phase symmetry is dynamically broken in the process of phase 
transition to the coherent regime. The role of the phason mode (p is to recover such a 
symmetry, thus re-establishing the gauge invariance of the theory. This is done through the 
emergence of the coherence implied by the phase locking between the matter field and the 
e.m. field. Let us see how this happens. 



B. The phase locking 

As shown above, provided 9q > ^, a. time-independent amplitude U{t) = U is compatible 
with the system dynamics (e.g. the ground state value of Aq in eq. (3.21) implies U = const., 
as noticed above). Eqs. (3.12) and (3.13) with am = a = di{t) — 5o(^) — ^{t) are 

U(t) = 2nAo(t)Ai(t) cos a(t) , (3.28) 
m = 2n^^^^§^sinait). (3.29) 

We thus see that U{t) = 0, i.e. a time-independent amplitude U — const, exists, if and only 
if the phase locking relation 



holds. In such a case. 



a^5,{t)-So{t)-ip{t)^^ (3.30) 



ip{t) ^ Si{t) - So{t) ^ u . (3.31) 



which shows that any change in time of the difference between the phases of the amplitudes 
ai{t) and ao{t) is compensated by the change of the phase of the e.m. field. When eq. (3.30) 
holds we also have Aq — — Ai (cf. eqs. (3.10), (3.11)). Provided ^0 > f > the phase relation 
(3.30) can be thus regarded as a further constant of motion implied by the dynamics: a = 0. 
It expresses nothing but the gauge invariance of the theory. Since 60 and (p are the NG 
modes, eqs. (3.30) and (3.31) also exhibit the coherent feature of the collective dynamical 
regime: the system of N dipoles and of the e.m. field is characterized by the "in phase" 
dynamics expressed by eq. (3.30) (phase locking). In other words, the gauge invariance of 
the theory is preserved by the dynamical emergence of the coherence between the matter 
field and the e.m. field. In such a regime we have 

aI-AI^ cos^ ^0 - - sin^ ^0 + 2C7V , (3.32) 
3 

to be compared with Alit) — Al{t) ^ at the thermal equilibrium in the absence of the 

collective dynamical regime discussed here. Eq. (3.32) shows the relevant role played by 

— 2 

the occurrence of a time-independent e.m. amplitude U ; the collective dynamical regime. 
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which sets in for Oq > j, allows that a non-zero permanent polarization appears when an 
electrical field is applied, as discussed in deriving eq. (2.14). In the following we will come 
back on this point. 

In conclusion we recognize that, starting ai t — from the initial condition |m(0)P = 0, 
and, correspondingly, from the zero order parameter uq — 0, a non-zero time-independent 
e.m. amplitude can develop (phase transition), provided 6q > as an effect of the radiative 
dipole-dipole interaction. This results in turn in the phase locking (3.30) and in the subse- 
quent coherence in the time behavior of the phase fields (cf. cq. (3.31)). Eqs. (3.30) and 
(3.31) show the role played by the phason field (f in recovering the gauge invariance in the 
process of phase transition to the collective dynamical regime. 

In the collective dynamical regime considered above the values of the amplitudes and 
Ai are related to the amplitude U through the relations (2.7) and (3.8). Moreover, we also 
obtain 

Al^^[l+ cos^ ^0 + (1 - I sin^ 2^o)'/'] ■ (3.33) 

which used in (3.18) shows that the oscillations around the ground state for Aq have pulsation 

i/ = 2^2(1 -isin2 2^o)'/^ 

The physical meaning of the phase locking can be stated as follows. The gauge arbi- 
trariness of the field A^ is meant to compensate exactly the arbitrariness of the phase of 
the matter field in the covariant derivative = — igA^. Should one of the two ar- 
bitrarinesses be removed by the dynamics, the invariance of the theory requires the other 
arbitrariness, too, must be simultaneously removed, namely the appearance of a well defined 
phase of the matter field implies that a specific gauge function must be selected. The above 
link between the phase of the matter field and the gauge of A^ is stated by the equation 

= dfj.^ (^/i is a pure gauge field). When (/^(x, i) is a regular (continuous differentiable) 
function then E = — ^ + VAq = ("^V + V^)(/9 = 0, since in such a case time derivative 
and the gradient operator can be interchanged. Analogously, in the space of the regular 
functions <^(x, t) it is B = V x A = V x = 0. Thus the existence of nonvanishing fields 
E and B in a coherent region implies that the time and space derivatives should act on a 
space larger than the space of regular functions, namely </7(x, t) should exhibit a (divergence 
or a topological) singularity within the region^. This is precisely what is observed, e.g., in 
type II superconductors when penetrated by the lines of a quantized flux in a vortex core. 

IV. DISCUSSION AND CONCLUSIVE REMARKS 

We have seen above that the rescaling of the field by the factor (cf. eq.(2.2)) induces 
the enhancement by the same factor of the coupling constants appearing in the field equations 
(2.15) (see also the coupling Vl introduced in eqs. (3.1)-(3.3)). This implies that for large 
the collective interaction time scale is much shorter (by the factor than the short range 
interactions among the atoms. Hence the mesoscopic/macroscopic stability of the system 
vs the quantum fiuctuations in the short range dynamics of the microscopic components. 
For the same reason, for sufficiently large N the collective interaction is protected against 
thermal fiuctuations. Indeed, thermal fiuctuations could affect the collective process only 
when kT is comparable or larger than the energy gap, whose value thus determines the 
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height of the protection; the larger the energy gap, more robust the protection. We do not 
present in this paper an estimate of the energy gap. We will do that in a future work. The 
role of the factor \fN in setting the time scale of the system can be understood also in 
connection with the influence on the system of atoms of an electric held E due, e.g., to an 
impurity (or to any other external agent). By closely following ref.^^, we will see indeed that 
for large the system of atoms behaves as a collective whole. The interaction H. — — d • E 
(eq. (2.12)) with the electrical field can be written^^ as 

H = h-i{h'^(T- + 6(7+) , (4.1) 

which is a Jaynes-Cummings-like Hamiltonian, as already mentioned in Section I in connec- 
tion with Haken analysis. In eq. (4.1) 7 is the coupling constant which is proportional to 
the matrix element of the atomic dipole moment and to the inverse of the volume square 
root h is the electric field quantum operator, are the atomic polarization op- 

erators. Let |0)j and |l)j, i — 1,...,N, denote the ground state and the excited state of 
each of the N two-level atoms, respectively, associated to the eigenvalues =f| of the op- 
erator (Jsi = |(|l)ii(l| — |0)ii(0|), no summation on i (see e.g. Section III. 6 of ref.^^ for 
the formal equivalence of the system of two- level atoms with a system of ^ spins). The 
operators (j+ = |1)m(0| and a~ — ((J+)^ generate the transitions between the two lev- 
els induced by the action of the electric field. The A^-atom system is thus described by 
(7+ = J2i^=i (^ti c^s = Sill c"3i with the fermion-like su(2) algebra 

[a3,(7±]=±a±, [(7-,a+] = -2(73. (4.2) 

Suppose that the electric field action induces the transition |0)j — > for a certain number 
of atoms, say / (as far as iV 3> / our conclusions will not be affected by the fact that initially 
some of the atoms are not in their ground state) . The system state may be then represented 
as the normalized superposition \l) given by 

10 = (|0)i|0)2...|0)Ar_i|l);v-m|l)7V-/+2...|l)7V + ■ ■ ■ 

+|l)i|l)2...|l),|0),+i|0),+2...|0)Ar)/y^. (4.3) 

The difference between the number of atoms in the excited state and the ones in the ground 
state is measured by (73: 

{l\a,\l)^l-\N (4.4) 

and the non-zero value of this quantity (proportional to the system polarization) signals 
that the dipole rotational {SU{2)) symmetry is broken. Operating with cr^ on |/) gives: 

(7+10 = ^/T^^^W^l\l + 1) , 

(7-|0 = ^N-{l-l)Vl . (4.5) 

Eqs. (4.4) and (4.5) show that (73 and are represented on |/) by 
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(73 = S+S 

^-^'^ " 'l-^S-, (4.6) 




N V N ' y/N \ N 



where = {S+Y , [S-,S+] = 1, S+\l) = VTTT\l + 1) and S-\l) ^ VJ \l - 1), for any 
I. Eqs. (4.6) are the Holstein-Primakoff non-hnear boson reahzation of 5'C/(2)^^'^^. ^ in 
eqs. (4.6) still satisfy the su(2) algebra (4.2). However, since ior N ^ I eqs. (4.5) give 

^ 10 = \l) , (4.7) 

the su(2) algebra (4.2) contracts in the large N hmit to the (projective) e(2) algebra (or 
Weyl-Heisenberg algebra) ^^■^'''■^^ 

[Ss,S^]=±S^, [S~,S+] = 1. (4.8) 

where 5*3 = a^. From eqs. (4.7) and (4.8) we see that, for large A^, denote the creation 
and annihilation boson operators associated to the quanta of collective dipolc waves excited 
by the electric field. The interaction (4.1) can be now written in terms of as 

H = hVNj{b^S- + bS+) , (4.9) 

We thus conclude that in the large N limit the collection of single two-level (fermion-like) 
atoms appears as a collective bosonic system. The original coupling of the individual atoms 
to the field gets enhanced by the factor ^/N and appears as the coupling of the collective 
modes S"^ (the system as a whole) to the field. We observe that, as shown by eq. (2.14), the 
polarization persists as far as r is non-zero, namely as far as the field E is active (i.e. H 0). 
The system finite size prevents indeed from having a persistent polarization surviving the 
— > limit^. In such a limit the dipole rotational symmetry is thus restored. 
Finally, we note that the collective dynamical features presented here are not substan- 
tially affected by energy losses from the system volume, which we have not considered in the 
discussion above. These losses are related with the different lifetimes of our different modes, 
according to the different time scales associated to the pulsations mo, cr and /i. An anal- 
ysis of energy losses when the system is enclosed in a cavity has been presented elsewhere 
in connection with the problem of efficient cooling of an ensemble of N atoms^^. Another 
problem which we have not considered in this paper is the one related to how much time the 
system demands to set up the collective regime. This problem, which is a central one in the 
domain formation in the Kibble-Zurek scenario, will be the object of our study in a future 
work. Here we remark only that, since the correlation among the elementary constituents is 
kept by a pure gauge field, the communication among them travels at the phase velocity of 
the gauge field. 
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